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Abstract— Many of us are interested in detecting the 
irregularity in very small region of a transient signal 
which cannot be detected by necked eyes. That can be 
possible by discretizing the region in which the 
irregularity lie. There are two ways to discretize the 
signal. In theory one can use Fourier tran s form on the 
signal and cut in M-pieces. Where as in practice 
discretization can be obtained by applying multiple 
number of filters (Filter Bank). 

This paper will give knowledge of wavelets and filter 
banks and later on connection between them. These 
are rapidly developing topics in real time. The 
technique of filter banks (for discrete signals) and 
wavelets (to represent functions) are used throughout 
signal and image processing for compression, 
denoising, enhancement, motion estimation and 
pattern recognition. New wavelets continue to be 
constructed for new applications. 

To understand this process we need to have a brief 
knowledge about the basic concepts like Fourier 
Transform, Short-Time Fourier Transform, Wavelet 
Transform, MRA, Filter and Filter bank. 

Fourier Transform: 

A Fourier transform of f (t) G L 1 (R) is given as 
F(w) = F(/(t)) = f^e~ iwt f(t)dt 
= <f(t), e iwt > 

It is also hold for f (t) G L 2 (R) 

Fourier transform is a mathematical description of the 
relationship between function of time and 
corresponding function of frequency. 

Fourier transform widely used in science and 
engineering that convert a signal which is intensity vs 


time in to a spectrum which is an amplitude vs 
frequency. 

Limitation of Fourier Transform: 

The following are limitations of Fourier Transform. 

1. It can be computed for only one frequency at a 
time. 

2. A small change in the signal would affect the 
entire frequency spectrum of a signal. 

3. The formula does not give any information about 
frequencies which develops with respect to time. 

Therefore, to overcome above limitations we switch 
on to Short-Time Fourier Transform. 

Short-Time Fourier Transform (STFT): 

A STFT is used to get approximate information from 
both time and frequency domain simultaneously. It is 
defined as 


STwf(to,w)= I f(t)w(t — tci)e lwt dt 

j — 00 

= <f (t), w (t-to) e lwt > 

Now we shall show how to discretize the signal by 
using STFT. 
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The STFT represents a sort of compromise between 
the time and frequency based views signal. It provides 
some information about both when and at what 
frequencies a signal event occurs. However, you can 
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only obtain this information with limited precision, 
and that is determined by the size of the window. 

Limitation of STFT: 

The limitation is that once you choose a particular 
size for the time window, that window is the same for 
all frequencies. Many signals require a more flexible 
approach one where we can vary the window size to 
determine more accurately either time or frequency, 
but that is possible in wavelet transform. 

Wavelet: 

History of wavelet is in 1983, Jean Morlet a French 
Geophysicist working in the oil company first 
introduce the concept of a wavelet. He developed a 
new time frequency analysis calling “Wavelets of 
constant shape.” 

In 1984, Grossman and Morlet said the mathematical 
foundation for the wavelet transform using the tools 
like Hilbert space, frame etc. 

In 1986 Daubechies, Grossman, Meyer developed the 
theory regarding wavelets. This was the beginning of 
the wavelet or rather golden era of development of 
wavelet concept. 

We all are interested in knowing what wavelet is. 
How it is useful. Basically small waves are called 
wavelets. 

Wavelet has small concentrated burst of finite energy 
in the time domain (this condition makes the wavelet 
little) and it exhibits some oscillation in time (this 
conditions makes wavelet wavy), hence the word 
wavelet. 

Thus, wavelets have oscillatory behavior in some 
interval and then decays rapidly to zero, outside the 
interval. 

One major advantage afforded by wavelets is the 
ability to perform local analysis that is, to analyze a 
localized area of a larger signal. 

Consider a sinusoidal signal with a small 
discontinuity one so tiny as to be barely visible. Such 
a signal easily could be generated in the real world, 
perhaps by a power fluctuation or a noisy switch. 


s ’ 0 Sinusoid with a small discontinuity 

' ' ; I 


A plot of the Fourier coefficients of this signal shows 
nothing particularly interesting; a flat spectrum with 
two peaks representing a signal frequency. However, 
a plot of wavelet coefficients clearly shows the exact 
location in time of the discontinuity. 
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Wavelet analysis is capable of revealing aspects of 
data that other signal analysis techniques aspects like 
trends, breakdown points, discontinuities in higher 
derivatives and self-similarity. Furthermore because it 
affords a different view of data than those presented 
by traditional techniques, wavelet analysis can often 
compress or de-noise a signal without appreciable 
degradation. 

Indeed, in their brief history within the signal 
processing field, wavelets have already proven 
themselves to be an indispensable addition to the 
analyst’s collection of tools. Hence it is very popular 
in field of analysis. 

Equate wavelets with sine waves, which are the basis 
of Fourier analysis. Sinusoids do not have limited 
duration they extend from minus to plus infinity. And 
where sinusoids are smooth and predictable, wavelets 
tend to be irregular and asymmetric. 

WA 

Sine Wave Wavelet (dblO) 

Fourier analysis consists of breaking up a signal into 
sine waves of various frequencies. Similarly, wavelet 
analysis is the breaking up of a signal into shifted and 
scaled versions of the original (or mother) wavelet. 
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As in above image of sine waves and wavelets, seen 
that intuitively that signals with sharp changes might 
be better analyzed with an irregular wavelet than with 
a smooth sinusoid, just as some foods are better 
handled with a fork than a spoon. 

Wavelet Transform: 

Wavelet transform represents a windowing technique 
with variable-sized regions. Wavelet transform allows 
the use of long time intervals where we want more 
precise low-frequency information and shorter regions 
where we want high frequency information. 

Here what this looks like in contrast with the time- 
based, frequency-based and STFT views of a signal: 






& 






o> 

=1 

C T 






m 

Ll_ 



Time Amplitude 

Time Domain (Shannon) Frequency Domain (Fourier) 


Time 

STFT (Gabor) 


£ 


Time 

Wavelet Analysis 


As in above image noticed that wavelet transform 
does not use a time-frequency region, but rather s 
time-scale region. 

The Continuous Wavelet Transform (CWT): 

If \|/ C L 2 (R) satisfies the “admissibility” condition 


C(p = (^) dt ; /C L 2 (R) 


= L 00 /( t )vb f a(t) dt 


= <f(t), V|/ b ,a (t) > 


where vpb.a (t) = 






Where a,bG R with af 0 as parameter a,b gives effect 
of magnification or contraction and translation 
respectively. \|/b, a (t) are called baby wavelets. 


Wavelets are very very small and having irregular 
behavior in very small interval. Therefore area 
covered under the curve is approximately zero. 

i.e. q£>(t)dt = 0; (p e L 2 (R) where 

(pit) is the basic or mother wavelet. 


Scaling: 

Wavelet analysis produces a time-scale view of a 
signal, and now what is exactly mean by scale. 


Scaling a wavelet simply means stretching (or 
compressing) it. To go beyond colloquial descriptions 
such as “stretching” that is introduce the scale factor, 
often denoted by the letter a . In sinusoids, the effect of 
the scale factor is very easy to see: 



The scale factor works exactly the same with 
wavelets. The smaller the scale factor, the more 
“compressed” the wavelet. 



fit) = y(f) ; a = I 

fit) = y(2f) ; a = | 

fit) = v(4f) ; a = 1 


It is clear from the diagrams that, for a sinusoidal sin 
(wt), the scale factor a is inversely propotional to the 
radian frequencyo>. Similarly, with wavelet analysis, 
the scale is inversely proportional to the frequency of 
the signal. 

Shifting: 

Shifting a wavelet simply means delaying its onset. 
Mathematically, delaying a function f(t) by k is 
represented by f(t-k). 

Discrete Wavelet Transform: (DWT) 

The continuous wavelet transform has unique 
advantage as its widely width is control by scale ‘a’. 
However computational load of continuous wavelet 
transform is very heavy in order to capture all 
characteristics. To minimize excess of computation in 
continuous wavelet transform mathematician has 
developed discrete wavelet transform which can be 
obtained from the process known as multi resolution 
analysis. 

Multi Resolution Analysis: (MRA) 

MRA decomposes space L 2 (R) in a set of 
approximate subspace {Vj} j 6 Z 

Vj =.® Wj_ 2 © Wj_ x 
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Where { Vj } is generated by some scaling function 
O 6 (I?) and {Wj } is generated by some wavelet 
i peL 2 (R) 

Which satisfies the following axioms: 

1. {0} c •••.T_ 2 c VLi c 

•••.c L 2 (R) therefore, the 

sequence of subspaces Vj is nested. 

2 . r\ Jez Vj = { 0 } 

3. Close U jezVj = L 2 (R)i.e. closure of 
union of Vj gives whole of L 2 (R) 

4. V J+1 = Vj © Wj 

5. fix) e Vj ^ /( 2x) e Vj+i, jez 

6. There exist a scaling function 
0(x)e Vq such that {<D ;fc } is an 
orthonormal basis that spans Vj. 

7. F(x) 6 T 0 /( 2x) e 

Vj+i, jez 

Filter: 

Filter is an operator which maps one signal to another 
signal. 

Filter has basically two properties. 

> Linearity: 

If H is a filter? Operator. 

H((a/+ fig) = aHif) + PH(g) 

Where 

a, f> are constants and f, g are two functions if, 

I) 

> Time invariant or Shift invariant: 

A filter H is called shift invariant or time invariant if 
it communicates with the shift operator HD = DH. 
This means that delaying the filtered signal is the 
same as filtering the delayed signal. 

■ When filter is operated on an impulse say S 
then the effect generated is known as impulse 
response. The impulse response is given by 
filter coefficients. 

h = (h n ) = H 5.(1) 

■ The fourier transform of the impulse response 
is known as frequency response of the filter. 

H(e lw ) = Ym hn e~ inw ( by using z — 
transform and taking z = e lw ) 

But frequency response is the combination of 
amplitude and phase angle so, 

H(e iw ) = e - ‘ nw |e i<l>(w) 


= |H(eiw) |e i<t>(w) 

We can describe filter in three different way that is in 
time domain, z-domain and in term of linear algebra. 

Filter Bank: 

Suppose discrete or an analog signal ‘s’ which is 
limited with bandwidth Q. We can write signal as the 
sum of M signals. Each of which have bandwidth 
Q/M. In this way, a wide band signal can be split into 
M signals of smaller band and transmitted over a 
channel with smaller bandwidth. The receiver can 
reconstruct the original signal. In theory we can 
compute, Fourier transform of signal ‘s’, cut this in M 
pieces and back transform. In practice, we are 
applying M filters to the signal and each of these 
filters generates one signal with limited bandwidth. 
This is called M channel filter bank. 

Now, suppose we are applying 5 filters to the signal 
and each of these filters generates one of the 5 signals 
with small bandwidth. This is called 5-channel filter 
bank. 

Now M=2 

Suppose a discrete signal 4 s’. We can apply a low pass 
filter H and a high pass filter G. which splits the 
signaTs’ in two parts. The part which contains the low 
frequencies which gives a low resolution idea of the 
signal and the other part which contains the high 
frequencies which gives the detail information. This is 
called 2-channel filter bank. It splits frequency band 
in two sub bands. 

It is possible to stretch the half bands again to the full 
bandwidth by down-sampling. Down sampling means 
signals passing through filter in which we delete odd 
samples and keep only even number of samples. 
Down sampling is denoted by (i) 

If s=s n is given signal. 

S’ = l s if s n ‘ = s 2n 

In general, s’ = (i M )s if s n ‘ =s„m (M is number of 
filters). 

In z-domain, 

S’= l s osV)=^ 

(by applying low pass filter & high pass filter G(z) = 
H(-z) ) 

„ v , + H i-z 1 / 2 ) 

S (Z) = --- 
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S (z ) 


5(z x / 2 ) + 5(—z 1 / 2 ) 
2 


In the frequency domain, 


S'(e iw ) 


5 ( e iw/ 2 ) + 5 ( e -i /2) 
2 


This shows that if bandwidth of S is 7i then band 
width of S’ is 2 ji. 


Thus the operation in upper branch of the analysis 
part gives LP (z) as the result of multiplication with 
the adjoint of the matrix. 


H = 


hS h3 
h 4 


hi h_ t 
h 2 K 
h 3 h x 


-(*) 


For the high pass band rc/2 < \w\< n we have to shift 
the spectrum first to low pass band 


|w|<jt/2 which corresponds to adding 7i to w. This 
means that for the shifted spectrum. The frequency 
response is given by 



(W + 7T) __ 


^s n e" £ e~ inn 

V (-1 ) n s n e- inw 


On the analysis side of a two channel filter bank. We 
have the application of the two filters H~* and G * 
which are both followed by down sample operation. 
Here we use the notation H * to denote that the 
transfer function of this filter is H *= £/< h k z k 


Similarly for G * to denote that transfer function of 
this filter is G *= £/c g k z k H * (z) is the transform of 
the time reversed sequence h where h = hk 

On the synthesis side, of the filter bank one finds the 
mirror image of the analysis side. 

First the signals are up sampled. This means that 
between every two samples a zero is introduced. This 
is denoted as 


s' = T s s 2n = s n and s 2n _ 1 = 0 S (z) 

= S(z 2 ) 

After the up sampling some filters H and G are 
applied and two result signals Y and X are added. The 
synthesis side should up to the analysis computing 
such that the resulting signal S’ is again equal to the 
original signal S. 


The operation of the filter bank can also be written in 
terms of (infinite) matrices. Filtering by low pass 
filter H * means multiplication with the infinite 
toeplitz matrix with entries the impulse response h* = 
(h - n ) i.e. the Fourier coefficient of H *. 


Down sampling means that we skip every other row 
(the odd ones). 


Similarly, for the lower branch. Thus the vector p of 
the samples of LP (z) and the vector q of the samples 
of HP (z) are obtained from the samples s of S(z) as 
follows 



On the synthesis side, the up sampling followed by 
filtering with H means that we multiply with toeplitz 
matrix whose entries are the impulse response 
coefficients (i.e. the Fourier coefficients of H(z)) and 
in which every other column is deleted. That is the 
matrix H which is defined like H but without the 
tildes. That is as follows 




The matrix G can be defined similarly for the other 
branch on the synthesis side. 

The samples s of the result 5 (z) are then computed 
from the samples of LP (z) and HP (z) by 

s = H p + G q 

= [h g; 



We shall have s = sifK*K = I 


The recursive application of a two channel filter bank 
leads to an M channel filter bank. 


If the 2-channel filter banks split in equal band-widths 
then the bandwidths of the end channels will not be 
the same. 


Wavelet Decomposition and Reconstruction: 

Consider the general structure of multiresolution 
analysis and wavelets as in 
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Vj =. ®Wj- 2 eWj- lt jez 

Where {Vj} is generated by scaling 
function® e L 2 (R) and {Wj } is generated by some 
wavelet ip e L 2 (/?). 

By third property of MRA we know every function f 
in L 2 (R) can be approximated closely as is desired by 
f N eV N , for some N e Z 

As, Vj = Vj-i© Wj_-L , for any j6 Z 

f N has a unique decomposition f N = f N -i © g^-i where 
/w-i 6 Vn-i an d 9 n -i 6 W N _i 
by repeating this process, we have 

/v = 9n -i + 9n-2 + f 9n-m + /n-m - 

-- ( 6 ) 

Where fjeVj and gjeWj and M is so chosen that f N - M 
is sufficiently “blurred”. The decomposition in (6) is 
unique is called “wavelet decomposition” ; and the 
“blur” is measured in terms of the variation( or 
frequency or number of cycles per unit length) of 

/w-M- 

A less efficient “stopping criterion” is to require 
H/n-mII to be smaller than some threshold. 

Now an algorithm approach for expressing f N as a 
direct sum of its components gN-i> 9 n- 2 > — > 9n-m 
and f N - M and recovering f N from these components 

Since both the scaling function Of Vo and the wavelet 
ip e W 0 are in V i and since V i is generated by 

® lfc (x) = 2zO(2x — k); k EZ 

There exists two sequences {pk} and {qk}£ l 2 such 


that 

®(x) = SfeP fc O(2x-k)-(7) 

<p(x) = S fc q fc O(2x-k) - (8) 

For all x G R 


The formulas (7) and (8) are called the “two scale 
relation “ of the scaling function and wavelet 
respectively. Since both 0(2x) and 0(2x — 1) are in 
Vi and Vi = V 0 © W 0 

There are four l 2 sequences which are denoted by {a. 

2 k}?{b- 2 k} 5 {^ 1 - 2 k} 5 {bl- 2 k}> k £ Z 

Such that 

0(2x) = X fc [a_ 2fc O (x - k) + b_ 2k <p(x - k)]- 

-(9) 


0(2x - 1) = X fc [ai_ 2fc O (x - k) + b^k <p( x - 

k/—- (10) 

This two formulas (9) and (10) can combine into a 
single formula: 

0(2x - 1) = (x - k) + b t _ 2k <p(x - 

k ],IEZ (11) 

Which is called ‘decomposition relation’ of O and (p. 

Now two pairs sequences ({pk}, {qk}) and ({ak}, {bk}) 
all of which are unique due to direct sum relationship. 
V! = Vo© W 0 

These sequences are used to formulate the following 
reconstruction and decomposition algorithm. Here 
{pk} and {qk} are called reconstruction sequences, 
while {ak} and {bk} are called decomposition 
sequences. 

CONCLUSION 

The lifting scheme which is an efficient 
computational scheme to compute wavelet transform 
can be derived directly from the polyphase matrix; it 
can be applied in much more general situation than 
classical wavelet fdter banks. Therefore, we introduce 
it via an alternative approach to wavelets (subdivision 
schemes) which will give a framework in which it is 
easier to consider more general situations. 

The signal is split into a low pass and a band pass part 
by the filters H and G respectively. This corresponds 
to computing the Vn and the Wn part of a signal in 
Vn+i. The results are then sub sampled, to remove 
redundancy. Next a (primal) lifting step is executed 
with the filter. At the synthesis side the same 
operations are undone in opposite order to obtain 
original signal. 
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